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C\ ! Abstract 
On 

' Quantum open systems are described in the Markovian limit by master equa- 

■ tions in Lindblad form. I argue that common "quantum jumps" techniques, 

. which solve the master equation by unraveling its evolution into stochastic 

£NJ ! trajectories in Hilbert space, correspond closely to a particular set of deco- 

| herent histories. This is illustrated by a simple model of a photon counting 

experiment. 

Recently a great deal of work has been done in quantum optics on "quantum jump" 
simulations of continuously measured systems with dissipation In this technique, a 

system described by a master equation for the reduced density operator p in the Markovian 
approximation ||, 

. ! 

p = -i[H, p\+J2 LmpLl - -L ] m L m p - -ptf m L mi (1) 



is "unraveled" into a jump process for pure states. H is the system Hamiltonian, and the 



{-L m } are a set of Lindblad operators which model the effects of the environment. 

Around the same time, the decoherent histories formulation of quantum mechanics was 
developed [|7|-[T2[]. In this formalism, one describes a quantum system in terms of an ex- 
^ ■ haustive set of possible histories, which must satisfy a decoherence or consistency criterion. 

Histories which satisfy this criterion have probabilities which obey the usual classical prob- 
ability sum rules. 

Both quantum trajectories and decoherent histories describe a quantum system in terms 
of alternative possible evolutions; they thus bear a certain resemblance to each other. What 
is more, quantum jumps are commonly interpreted as giving the results of continuous mea- 
surements; and histories which correspond to records of a "classical" measuring device should 



always decohere |10[. Thus, there should be a set of decoherent histories which correspond 
to the quantum trajectories of a continuously measured system. 

Exactly such a correspondence has recently been shown between decoherent histories and 
quantum state diffusion (QSD), another unraveling of the master equation, by Diosi, Gisin, 
Halliwell and Percival Though this result was pioneering, it was rather abstract, and 
lacked any direct connection to a physical measurement situation. Similar results for yet 
another unraveling were given by Paz and Zurek and Diosi [[Tgl , in a model with exact 
decoherence, but also far removed from physical measurement situations. Other treatments 
T6| have been framed in terms of measurement alone. 



1 



Consider a quantum system with a Hamiltonian H , completely isolated except for a 
single channel of decay, which is monitored by an external photon detector. We model this 
detector as a single two-level system (the "output mode") with states |0) and |1) strongly 
coupled to an environment representing the remaining degrees of freedom of the device. 

The measuring device produces two important effects. The first is dissipation. Excita- 
tions of the output mode will be absorbed by the measuring device with a rate Ti which 
we assume to be rapid compared to the dynamical time-scale of the system. The time 1/ri 
represents the time-resolution of the detector. 

The second effect is more subtle but just as important: decoherence. As the state of 
the output mode becomes correlated with the internal degrees of freedom of the measuring 
device, the phase coherence between the ground and excited states of the output mode is 
lost. Investigations of this process have shown that the loss of coherence is generally far 
quicker than the actual rate of energy loss |T1J. This decoherence rate is T 2 IY 

We suppose that the system is linearly coupled to the output mode via the Hamiltonian 

Hi = k(ci ] <g>6 + a<g>& f ), (2) 

and the total Hamiltonian is 

H = H ®i + K(a^ ®b + a®ti) } (3) 

where a and b (a 1 and o 1 ) are the lowering (raising) operators for the system and output 
mode, respectively. The hierarchy of evolution rates is T 2 ^> T\ ^> k. 
The total system obeys the master equation 

p= -i[H, p\ + FtbpP - Y^ b P ~ Y P ^ b 
+T 2 a z pa z - T 2 p = Cp, (4) 

where p is the density matrix for the combined system and output mode, and the Pauli 
operator a z acts on the output mode. £ is the Liouville superoperator. This is a linear 
equation, and so can be formally solved: 

pfo) = exp{cfo-t 1 )\p(t 1 ). (5) 

Assume that we start in a pure state \^/) = \ip) (g> |0). We can expand p 

p(t) = PooOO ® |0> <0| +p 01 (t) (8) \0)(1\ + p w (t) ® |1)(0| +p n (t) ® |1)(1|. (6) 

In terms of these components the master equation becomes 

poo = -i[H , Poo] - iKalpio + inp 01 a + Tipn, 
Poi = -i[H , Poi] - ina'pn + inp^a) - Gp m = p\ , 

p n = -i[H , p n ] - map i + iKp 10 a^ - Tipn, (7) 

where G = r\/2 + 2T 2 T± ^> n. (This combination G occurs frequently in the equations 
which follow.) 
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Since the poi,Pio,Pn components are heavily damped, we can adiabatically eliminate all 



components other than p 00 |LJ 



2k 2 k 2 , K 2 

Poo = -i[H , Poo] + -^-apoo« f - -g^ap 00 - —p 00 a i a, (8) 

to first order in k 2 /G, provided that the system is not so highly excited as to emit too 
rapidly, i.e., /t(a^a) <C IY 

We can unravel the master equation (f|) into a sum over quantum jump trajectories. 
First, define a non-Hermitian effective Hamiltonian 

H cS = H - i(K 2 /G)a)a. (9) 

Assume that the system (excluding the output mode) begins in a pure state \ip). \ip) evolves 
according to the Schrodinger equation, 

^ = -i*-i*>. do) 

interrupted at random times by sudden quantum jumps 

M-aM- (11) 

These jumps correspond to the detection of photons Note that this evolution does not 
preserve the norm of the state. The physical state is taken to be = / <J {ip\ip), the 
renormalized state. 

The probability that an initial state evolves for a time T and undergoes N jumps 
during intervals 5t centered at times t±, . . . , tjv is 

(2StK 2 /G) N Ti^e~ iA ^^ T ~ tN ^ae- lfl ^ tN - tN ^a ■ ■ ■ ae"^* 1 

x l^le^Wt • • -aV^( T -^)), (12) 



i.e., the norm of the unrenormalized state gives the probability for that state to be realized. 

Equation (|Sp is valid only as long as the Markovian approximation remains good. In 
the case of our toy model, this means that it is valid only on time-scales longer than l/IY 
Thus, rather than a jump occurring at a time tj, it is more correct to consider the jump as 
occurring during an interval 5t ~ 1/ri centered on tj. This is fine so long as the jumps are 
separated by more than 5t on average, i.e., the system is not too highly excited. 

By averaging \ip)(tp\ over all possible trajectories with the probability measure (0), one 
can show that this unraveling reproduces the master equation (||) as required 0. 

Now, let us turn to the decoherent histories picture. In non-relativistic quantum mechan- 
ics, a set of histories for a system can be specified by choosing a complete set of projections 
{V^.(tj)} at a sequence of times t±, . . . ,t^, which represent different exclusive alternatives: 

E K (*i) = i> K, (^)K' (fs) = ( 13 ) 
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A particular history (denoted h) is given by choosing one V at each point in time. The 
decoherence functional on a pair of histories h and h! is 

D[h, h>] = Tr[vl(t N ) ■ ■■Vl 1 {t l )p{t )V l a ,{t l ) ■ ■ -P^}, (14) 



where p(to) is the initial density matrix of the system ||10|| . This satisfies the decoherence 
criterion if the off-diagonal terms vanish, D[h, h'] — 0, h ^ h! . The diagonal terms then give 
the probabilities of the histories, p(h) = D[h, h]. 

Suppose our initial pure state is \^) = \ip ) <S> |0), and we consider histories composed 
only of the Schrodinger projections 

A, = i®|o>(o|, A = i®|i><i|, (is) 

representing the absence or presence of a photon in the output mode. The projections are 
spaced a short time 5t apart, and a history is composed of N projections, representing a 
total time T = N5t. A single history h is specified by a string {oti, a 2 , ■ ■ ■ , a^}, where 
a,- = 0, 1. In this case, by the quantum regression theorem EtJ the decoherence functional 



(|14]) becomes 

D[h, ti] = Tr|p ajv e £5i (P Qjv _ ie £5i (- • • e m (V ai |*> ^\P a[ ) ■ ■ .)^}. (16) 

The Liouville time evolution superoperators (|]) evolve pure states into mixed states. 
This is counteracted by the effect of the repeated projections P a . 

From the equations (^) we can determine the character of the different histories. The 
crucial parameter is the size of the spacing 5t between projections. The interesting regime 
is in the range 

^ « St « ^ (17) 

On this time-scale, the T 2 terms are sufficient to insure decoherence while the effects of the 
Ti terms are resolved into individual pure state trajectories. This last is a subtle point. The 
probability of a photon being emitted in any single time step is small. However, if a photon 
is emitted, it has an appreciable possibility of being absorbed on a time scale l/TY The 
effect of decoherence produces the terms (k 2 /G)d^ap 00 and {k 2 / G) p^a) d in equation (||), 
which are included in the effective Hamiltonian (^). These terms are already important on 
a time scale 5t l/IY By contrast, the term (2k 2 / G)dp 00 d^ is produced by the effects of 
dissipation, which only become important on a time scale l/IY It is this term which causes 
pure states to evolve into mixed states in equation (R). By choosing a time 6t l/IY we 
can maintain the purity of the system state over a full trajectory, as we shall see. 

If the external mode is initially unexcited, with p = p 00 eg) |0)(0|, then after evolving for 
a time St the state becomes 

K 2 K 2 

(e CSt p)oo= Poo - i[H , poo}St - —&ap 00 St - —p m d ] dSt + h.o.t. 



-i{H -iU 2 /G)a< a)St „ i(H +i(K 2 /G)^ a)8t 

e Poo6 j 



(^01=^00^ + h.o.t. = (e^p) f io, 
2k 2 

(e C5t p) u = —dpood^St + h.o.t. (18) 
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Here we see the appearance of the effective Hamiltonian H e s, just as in the quantum jump 
unraveling. 

If the initial state is p = pn <E> |1)(1|, after a time St the state becomes 

(e CSt p)oo= T^te-^pne^ + ^ p n a + h.o.t. 
(e £5 V)oi=-^tp 11 + h.o.t. = (e^p) t io, 

(e £<5 V)n= (1 - (ri + 2K 2 /G)5t)e- i ^ B&t p ll e iil l« &t + h.o.t., (19) 

Once again the effective Hamiltonian appears, together with two additional effects. The first 
is the possibility that the photon in the excited mode will be absorbed by the measuring 
device. The second (much smaller) effect is the possibility that the photon will be coherently 
re-absorbed by the system. This last process is so weak as to be negligible. 

By combining the above expressions with the appropriate projections Vq and V\ (which 
pick out the p o or Pn component, respectively), we can write down the probabilities of all 
possible histories. 

Note that the magnitude of the off-diagonal poi.io terms in both cases is of order 0(k/G). 
(This is also true for transitions from off-diagonal to diagonal terms.) This will be important 
in estimating the decoherence of this set of histories. 

First consider the history given by an unbroken string of N Vq projections, corresponding 
to no photon being emitted during a time NSt. 

The probability of such a history is given by the diagonal element of ( |TE| ) . We can expand 
the time evolution superoperator using ( |18|) and see that after St we get 

V e CSt (\il>)(ip\ ® |0)(0|)P « ^ e - i( ^°- i(K2/G)ata)5t |V')(V'|e 4( ^ o+i(K2/G)ata)<5t )|0)(0|. (20) 

Repeating this N times and taking the trace we get 

p(h) « Trje-^^IV)^!^^ 5 *}, (21) 

which exactly agrees with the probability of the quantum jump trajectory when no jumps 
are detected. 

Suppose now that at time NSt a photon is emitted, so that instead of using a final 
projection Vo we use V\. This corresponds to keeping the p n component of exp(CSt)p 
instead of p o, and yields a probability 

p(h) « (2StKyG)Ti^ae- i ^ mt \^}{ij\e i ^ mt a^, (22) 

Once again, this exactly agrees with the probability of the corresponding quantum jump 
trajectory. What happens after the output mode has "registered" as being in the excited 
state? Essentially, there are two possibilities: either the output mode can drop back to the 
unexcited state (representing absorption of the photon by the measuring device) or it can 
remain in the excited state. 

V e cst mm ® |1)(1|)P « rVtlVWI ® |0)(0|, (23) 
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V x e cst {\^){^\ ® « (1 - r 1 5t)e-^ fr<5 *|^ / )(V , '|e i ^ff <5t <g> |1)(1|. (24) 

We see that the output mode has a probability of roughly T^t per time St of dropping 
back to the ground state, while the system state continues to evolve according to the effective 
Hamiltonian H e s- 

This is slightly different from quantum jumps. Quantum jumps are resolved only on a 
time-scale l/r 1; not St -C l/IY However, there is a near- unity probability of the external 
mode returning to the ground state within a time of order l/r 1; so one can simply sum over 
all histories in which the photon is absorbed within this time. It is easy to see that these 
will, once again, match the quantum jump trajectories exactly. This type of coarse-graining 
is common in decoherent histories [pT0|, QJ] , and does not alter the form of the result. 



By combining the three cases described in this section, one can produce histories of 
multiple jumps. It is clear that the probability of such a history will be exactly of the form 
(0). 

In order for this discussion of probabilities to be meaningful we must require the histories 
to be decoherent. Exact decoherence is a very difficult criterion to meet. It is more usual 
to show that a model is approximately decoherent. In order for the probability sum rules to 



be satisfied to a precision e<lwe require that fll9 | 



\D[h,ti]\ 2 < e 2 D[h,h]D[h',ti] = e 2 p(h)p(ti), (25) 

for all distinct histories h, h! . Generally speaking, the more "different" a pair of histories is 
(i.e., the more projections they differ in), the more suppressed the off-diagonal term. So it 
suffices to look at two histories which differ at a single time tj, one having a projection Vo, 
the other V%. This is equivalent to picking out the poi or p w component of exp(CSt)\ip') 
at that time. 

Examining the components given by (plj|-|19|) , 



\D\h,h'}\ 2 1 , . 

1 (26) 



p(h)p(h>) (GSt) 2 ' 

we expect the sum rules to be obeyed with a precision of roughly 0(1/ GSt) (where we once 
again have assumed K(a^a) is small compared with Ti). 

We have seen how, in this simple model of a continuous measurement, the set of quantum 
jump trajectories corresponds to a set of decoherent histories. One of the principal goals 
of the decoherent histories program was to create a formalism which would reproduce the 
results of the usual Copenhagen formalism in measurement situations. It is pleasant to note 
that extensions to repeated or continuous measurements follow naturally within decoherent 
histories. 

In this letter, I considered only one measurement scheme: direct photodetection. In 
fact, there are many different schemes which give rise to different unravelings of the same 
master equation — heterodyne and homodyne detection, to name two @J21| . I have no doubt 



that arguments similar to those I have advanced in this paper will demonstrate similar 
correspondences to different sets of decoherent histories. 

This correspondence also has obvious practical benefits. Enumerating a full set of deco- 
herent histories and calculating their probabilities is an arduous and unrewarding task, in 
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general. There is a great deal of accumulated experience in simulating quantum trajectories; 
in situations where one would like to generate individual decoherent histories with correct 
probabilities, existing numerical techniques could be used. 

The decoherent histories formalism was developed largely in response to the problems of 
quantum cosmology, while quantum trajectories arose from problems in quantum optics and 
atomic physics. Both extend the von Neumann description of quantum mechanics to new 
realms of application. As the connections between the two formalisms are further explored, 
we can hope that a great deal of interesting physics will emerge. 

I would like to thank Lajos Diosi, Murray Gell-Mann, Nicolas Gisin, Jonathan Halliwell, 
Jim Hartle, Peter Knight, Mike Nielsen, Ian Percival, Martin Plenio, and Riidiger Schack 
for valuable conversations, suggestions, and feedback. After completion of this research, I 
became aware of related work by Ting Yu, from a rather different approach |22|. Financial 
support was provided by the UK EPSRC. 
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